Sudden death and sudden birth of entanglement in common structured reservoirs 
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We study the exact entanglement dynamics of two qubits in a common structured reservoir. We 
demonstrate that, for certain classes of entangled states, entanglement sudden death occurs, while 
for certain initially factorized states, entanglement sudden birth takes place. The backaction of the 
non-Markovian reservoir is responsible for revivals of entanglement after sudden death has occurred, 
and also for periods of disentanglement following entanglement sudden birth. 
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Entanglement is one of the most intriguing features of 
quantum mechanics [l| : hence understanding its proper- 
ties and dynamics is of paramount relevance for a num- 
ber of applications in modern physics. Furthermore, 
quantum cryptography, quantum information processing 
and quantum measurement are examples of branches of 
physics in which entanglement plays an essential role Q • 
Quantum properties, however, are very fragile: realistic 
quantum systems are not closed, and due to the interac- 
tion with the environment, their entanglement and coher- 
ence can be irretrievably lost. On the other hand, recent 
works have shown that entanglement can actually 

revive, or be preserved, using the quantum Zeno effect, or 
it can even be trapped. Such effects arise when the sys- 
tem of interest interacts with non-Markovian reservoirs. 
These reservoirs are characterized by structured spectral 
distributions [f|. 

In this Letter we propose an exactly solvable model for 
the time evolution of two qubits interacting with a com- 
mon structured reservoir. We investigate the entangle- 
ment dynamics for initially entangled states of the qubits, 
and also for initially factorized states. We identify the 
effects of the non-Markovian backaction of the reservoir 
during the disentangling process and in the entanglement 
birth phenomenon. In particular we prove that, due to 
the reservoir memory effect, entanglement sudden death 
(ESD) is enhanced as compared to the Markovian case. 
Contrary to a previous study on ESD in common envi- 
ronments our approach is exact and it does not rely 
on either the Born or the Markov approximation. 

The entanglement dynamics of two qubits interact- 
ing with two independent Lorentzian reservoirs are 
known [3| ■ In this Letter we focus on the common reser- 
voir scenario and compare it with the independent reser- 
voirs case. We note that, already in the Markovian case, 
the entanglement dynamics in common and independent 
reservoirs present striking differences. While two qubits 
interacting with two independent reservoirs can disen- 
tangle completely and permanently in a finite time @, Q , 
in a common reservoir entanglement can disappear for 
a finite time and then reappear again [loj . This is due 
to the fact that common reservoirs tend to create en- 



tanglement rather than destroy it completely, since they 
indirectly couple the two qubits [111 ]. As we show in the 
following, even when the qubits are initially prepared in 
a factorized state, the correlation created by the environ- 
ment can lead to a phenomenon which is known as the 
entanglement sudden birth (ESB) [lit ]. 

The analytic solution presented in this Letter is valid 
for a general initial state of the two-qubit system, and in- 
cludes states with two excitations. The exact solution for 
the single excitation case has been discussed in Ref. jij], 
where it was shown that the initial entanglement never 
disappears for finite periods of time, i.e., ESD never oc- 
curs. When two excitations are present, as for the model 
discussed here, the derivation of the exact analytical so- 
lution is much more complicated. We have proved, how- 
ever, that by using the pseudomode approach [131 ] and 
establishing a connection with a three-level ladder sys- 



tem [14J , it is still possible to solve the dynamics without 
performing any approximations. 

We consider a two-qubit system interacting with a 
common zero-temperature bosonic reservoir. Our cho- 
sen specific system consists of two two-level atoms inter- 
acting with the electromagnetic field. The initial state 
of the multimode field is the vacuum state. The Hamil- 
tonian of such a system in the rotating wave approxi- 
mation is given by H = Hq + Hi nt , which, in the basis 
{|00),|10),|01),|11)}, reads 
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Here, o~± and a± are, respectively, the Pauli raising and 
lowering operators for the atoms A and B, u>q is the Bohr 



frequency of the two atoms, a k and a\ are the annihi- 
lation and creation operators for the field mode k, and 
mode k is characterized by the frequency u>k and the cou- 
pling constant g^. For the sake of simplicity, in the fol- 
lowing we assume that the two atoms interact resonantly 
with a Lorentzian structured reservoir, such as, e.g., the 
electromagnetic field inside a lossy resonator. 
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Since the atoms are identical and equally coupled to 
the reservoir, the dynamics of the two qubits can be ef- 
fectively described by a four state system in which three 
states are coupled to the vacuum in a ladder configu- 
ration, and one state is completely decoupled from the 
other states and from the field. This can be shown by 
writing the Hamiltonian in the basis {|0) = |00), |+) = 
(|10) + |01»/V2, |-> = (|10) - |01))/V2, |2) = |11)} so 
that 



H = 2u Q \2){2\+u Q {\- 



+ |->H)+I>*4«*> ( 3 ) 



H int = ^2V2g k a k (\+){0\ + \2){- 



(4) 



The |+) and |— ) states are, respectively, the superradi- 
ant and subradiant states. From Eqs. and ^ it is 
apparent that the subradiant state does not decay. On 
the contrary, the superradiant state is coupled to states 
|0) and |2) via the electromagnetic field. 

The total Hamiltonian, given by Eqs. |3j) and ((U, ac- 
tually consists of two parts, a part describing the free 
dynamics of the |— ) state, and the remaining terms de- 
scribing a three-state ladder system ||0), |+), |2)} with 
the transitions |0) «-> |+) and |+) <-» |2) having the 
same frequencies and identically coupled with the com- 
mon bosonic reservoir. 

It is a mathematically and computationally demand- 
ing task to solve numerically the infinite set of differ- 
ential equations for the complex amplitudes appearing 
in the state vector of the total system. However, hav- 
ing m mind Eqs. © and (|4|), we can greatly simplify 
the analytical treatment by noting that the dynamics 
of the three-level ladder system, characterized by tran- 
sitions with equal frequencies and identically coupled to 
the same Lorentzian structured reservoir, can be exactly 
solved using the pseudomode approach [131]. In Ref. [15| 
this approach is extended to multiple excitations. 

The pseudomodes are auxiliary variables defined from 
the properties of the spectral distribution, and they al- 
low us to derive a Markovian master equation for the 
extended system comprised of the system of interest and 
the pseudomodes. Such an exact master equation de- 
scribes the coherent interaction between the system and 
the pseudomodes, and the latter, in turn, leak into inde- 
pendent Markovian reservoirs. 

The exact dynamics of the three-state system interact- 
ing with a Lorentzian structured reservoir is contained in 
the following pseudomode master equation 

dp r 

— = -i[V, p] - -{a)ap + pa) a - 2apa t ), (5) 

where p is the density matrix of the three - level system 
and the pseudomode, and 

V = V2n(o|+)<0| + a t |0)(+| + a|2)(+| + a f |+)(2|). (6) 



Here, a and a\ T and Vl are, respectively, the annihila- 
tion and creation operators, the pseudomode decay rate 
into its Markovian reservoir, and the coupling constant of 
the pseudomode to the ladder system. The pseudomode 
is associated with the Lorentzian spectral distribution 



JH = 



n 2 



r 



7T (W - LUo) 2 + (r/2)2 



(7) 



where F/2 describes the frequency width of the spectrum 
and is related to the reservoir correlation time. For small 
values of T the pseudomode can be associated with the 
real cavity mode of frequency u)q. 

In order to find the dynamics of the three states we 
solve the master equation in Eqs. (JSJ) and ([6]) and then 
trace out the degrees of freedom of the pseudomode. 
Then we take into account the |— ) state and consider 
the density matrix for the two atoms in the initial basis 
{|00),|10),|01),|11)}. 

Our aim is to investigate the effects of the non- 
Markovianity of the reservoir on atomic entanglement dy- 
namics. To quantify the entanglement we use the Woot- 
ters concurrence [lfjj], defined as C(t) = max{0, \f\[ — 
\/A2 — \f\l — V^h}, where {A,} are the eigenvalues of the 
matrix R = p((Jy®cr y 3 )p*(a y l (g)a y 3 ), with p* denoting the 

complex conjugate of p and <J^ B are the Pauli matrices 
for atoms A and B. This quantity attains its maximum 
value 1 for maximally entangled states and vanishes for 
separable states. 

For initial states of the form 



|tf(0)> = a|00) + e* 9 (l-a 2 ) 1/2 |ll) 



(8) 



the density matrix of the atomic system has an "X" form: 



P = 



( a(t) w(t) \ 

b(t) z(t) 

z*(t) c(i) 

\w*[t) d(t) J 



(9) 



with non-zero elements only along the main diagonal and 
antidiagonal. Due to the structure of the differential 
equations for the density matrix elements [see Eq. 
the "X" form is preserved during the evolution. Then 
the concurrence has a simple analytic expression 



C(t) = max{0,Ci(i),O>(i)}, 



where 



C 1 (t) = 2\w(t)\-2^b(t)c(t), 



C 3 (t) = 2\z(t)\ - 2^a(t)d(t). 



(10) 



(11) 



We plot the time evolution of concurrence as a function 
of both the parameter a 2 and the dimcnsionlcss quantity 
jot, where the parameter 70 = 4f2 2 /r is the Markovian 
decay rate of the atoms, i.e., the inverse of the atomic 
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Non-Markovian Common Rkskrvoir 




FIG. 1: (Color online) (a) Concurrence as a function of scaled 
time and a 2 for two atoms prepared in the state ((8| and inter- 
acting with a common Lorentzian structured reservoir. For 
comparison we show also the equivalent Markovian case (b) 
and the non-Markovian independent reservoirs situation (c). 



relaxation time in the Markovian limit [6(. The parame- 
ters T and f2 are also expressed in units of 70. We look at 
the dynamics in the strong coupling regime, which is ob- 
tained when r/fi < 4; in particular we choose T = 0.2 70 
and £1 = V0.05 70, corresponding to T/Q = y/0.8. These 
values are achievable experimentally, e.g., with circuit 
QED setups [nj]. 

Figure [TJ a) shows the behavior of concurrence for two 
qubits prepared in the state (jHJ ■ Depending on a 2 differ- 
ent dynamical behavior is clearly visible. For a 2 > 1/4 
the entanglement dynamics presents damped oscillations. 
For a 2 < 1/4 finite periods of complete disentanglement 
are followed by entanglement revivals. Entanglement re- 
vivals are amplified for stronger non-Markovian condi- 
tions, i.e., for smaller values of the ratio T/Cl. We have 
verified that when the ratio r/Q is, e.g., ten times smaller 
than the one used in Fig. [TJ the oscillations become much 
stronger in amplitude, the number of ESD periods and 
revivals increases, and they last for a longer time. On 
the other hand, when we increase r/fi, the ESD regions 
shrink and the period of the entanglement oscillations 
becomes smaller. We note in passing that entanglement 
oscillations also occur when the qubits share just one ex- 
citation (as shown in Ref. Q). In that case the memory 
depth of the reservoir leads to entanglement oscillations, 
but ESD does not occur. 

Entanglement sudden death and entanglement re- 
vivals, in the common structured reservoir, are basi- 



cally due to two combined and intertwined effects: the 
backaction of the structured reservoir and the reservoir- 
mediated interaction between the qubits. In order to un- 
derstand the role played by each of these effects we com- 
pare our results with the cases of the common Markovian 
reservoir, and two independent non-Markovian reser- 
voirs. 

In Fig. [UJb) we show the entanglement dynamics for 
two qubits in a common Markovian reservoir. As ex- 
plained in Ref. [lol ] , a period of complete disentanglement 
is followed by a revival of entanglement, but no oscilla- 
tions are present. Such a revival is due to the action of the 
common reservoir which tends to create quantum correla- 
tions between the qubits, providing an effective coupling 
between them. A comparison between Figs.[lja) andQJb) 
shows that the feedback of information from the reservoir 
into the system, characterizing non-Markovian dynamics, 
enhances the appearance of ESD regions, since it tends to 
recreate the conditions that led to the first ESD period. 

Figure [2c) shows the entanglement dynamics when 
the qubits interact with two independent non-Markovian 
reservoirs. Three different dynamical regions are clearly 
identified. Depending on the initial state of the system, 
entanglement can die after a finite time, or oscillate while 
going asymptotically to zero, or reappear after an ESD 
period. The revival phenomena in this case stem from the 
non-Markovian behavior of each single qubit interacting 
with its own reservoir [3}. No revivals of entanglement 
are present in the Markovian case § . 

A comparison between the non-Markovian cases of 
Fig. QJa) and QJc) reveals that, for the same type of 
reservoir spectrum, ESD regions are much wider in the 
independent reservoirs case than in the common reservoir 
case. Since both cases take into account memory effects, 
this suggests that the reservoir-mediated interaction be- 
tween the qubits, in the common reservoir scenario, effec- 
tively counters the fast disappearance of entanglement. 

Non-Markovian effects also influence strongly the dy- 
namics for initially factorized states, when the reservoir- 
mediated interaction between the qubits leads to entan- 
glement generation. As an example we take the initial 
state 

p(0) = (a 2 \0 A )(0 A \ + (l-a 2 )\l A )(l A \) 

®(a 2 \0 B )(0 B \ + (l-a 2 )\l B )(l B \). 

In the basis {|0), |+), |— ), |2)} this state takes the form 

p(0) = a 4 |0)(0| + (l- a 2 ) 2 |2)<2|) 

+ a 2 (l-a 2 )(|+)(+| + |->(-|). 

Again the density matrix has the "X" form so the con- 
currence can be calculated from Eqs. (Ti"0|) and (fTTj) . We 
evaluate the concurrence for the same parameters of the 
spectral distribution used in Fig.QJa). 

Figure [3 shows the evolution of concurrence when the 
atoms are prepared in the factorized state (fT2"]) . We note 
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FIG. 2: (Color online) (a) Concurrence as a function of time 
and a 2 for two atoms prepared in the factorized state (|12p and 
interacting with a common Lorentzian structured reservoir. 
For comparison we show the short-time region of (a) in (b), 
and the Markovian common reservoir result in (c). 



first that, for independent reservoirs, sudden birth of en- 
tanglement does not appear, while in our common reser- 
voir case, ESB does occur, as clearly demonstrated in 
Fig. [2^a). Moreover, in a structured reservoir such a 
phenomenon presents new interesting features, compared 
to the Markovian shown in Figs, db)-(c). To 

elaborate: in a Markovian reservoir, entanglement sud- 
den birth takes place at different times depending on the 
parameter a 2 . The smaller the value of a 2 , the longer 
the time taken for entanglement generation. In other 
words, the ESB time monotonically decreases with a 2 . 
However, in the non-Markovian case, ESB revivals occur 
and complicate the simple Markovian picture. In gen- 
eral, the reservoir memory prolongs the initial disentan- 
glement. In this case, therefore, the reservoir backaction 
dominates over the tendency of the common reservoir to 
create entanglement between the qubits. ESB periods 
and disentanglement revivals become more frequent and 
numerous for stronger non-Markovian conditions. 

It is worth noting that a necessary condition for ESB 
is that (— \p\— ) 7^ 0. Indeed, the long-time asymptotic 
value of concurrence is directly related to the subradiant 
state component of the initial state 12 1. 

Whatever the value of a is, for t — > oo the initial pop- 
ulation of the three state ladder system will decay to the 
|0) state, while the population of the subradiant |— ) state 
will be trapped. The asymptotic stationary state of the 
system has the form 



p(t-oo) = (l-fc)|0>(0|+fc|->H, 



(14) 



with k = a 2 (l — a 2 ). For a 2 ^ 0, or 1, this state is not a 
factorized state. The stationary value of the concurrence, 
calculated from the analytic solution, is C(t — ► oo) = k. 

In conclusion, we have presented a non-Markovian 
model describing the exact entanglement dynamics of two 
qubits interacting with a common structured reservoir. 
We have brought to light new entanglement features for 



qubits prepared in both entangled and factorized states. 
The non-Markovian nature of the reservoir protracts the 
disentanglement process while enriching the revivals, and 
at the same time it enhances the regions of ESD. The 
backaction of the reservoir slows down the generation of 
entanglement and further manifests itself in the appear- 
ance of periods of death and resurrection. The reservoir- 
mediated interaction between the qubits strikingly dis- 
tinguishes the dynamics in a common reservoir from the 
independent reservoirs case. Our predictions apply to 
cavity QED experiments with trapped ions, and to circuit 
QED experiments. In the first context, entanglement be- 
tween two remotely located trapped atomic ions has been 
recently demonstrated [l8| and multiparticle-entangled 
states can be generated and fully characterized via state 
tomography [191 ]. In the second context, field coupling 
and coherent quantum state storage between two Joseph- 
son phase qubits has been achieved through a microwave 
cavity on chip 17 , 0] ■ Due to the possibilities for re- 



alizing strong coupling conditions between atoms and a 
high finesse cavity [2l|, l22j a deep understanding of the 
non-Markovian dynamics is now indispensable. 

The authors thank CIMO, the Academy of Finland 
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